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Abstract 

The reparametrization transformation between ultrametrically organised states of repli- 
cated disordered systems is explicitly denned. The invariance of the longitudinal free energy 
under this transformation, i.e. reparametrization invariance, is shown to be a direct con- 
sequence of the higher level symmetry of replica equivalence. The double limit of infinite 
step replica symmetry breaking and n — > is needed to derive this continuous gauge-like 
symmetry from the discrete permutation invariance of the n replicas. Goldstone's theorem 
and Ward identities can be deduced from the disappearence of the second (and higher order) 
variation of the longitudinal free energy. We recall also how these and other exact state- 
ments follow from permutation symmetry after introducing the concept of "infinitesimal" 
permutations. 

1 Introduction 

Exact statements such as Goldstone's theorem [[j], ^] and Ward identities || proved very useful in 
the perturbative analysis of ordinary statistical systems with a continuous symmetry. For exam- 
ple, Goldstone's theorem ensures that the low temperature phase of the O(m) model is massless 
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(the transverse susceptibility is infinite), while the Ward identities, i.e. the exact relations be- 
tween the different vertex functions (derivatives of the Legendre-transformed free energy), help 
us to get rid of dangerous infrared divergences, order by order in perturbation theory Q. 

Zero modes and massless phases are surprisingly frequent in quenched random systems, too. 
The idea of the spin glass phase being marginal occured very early to various authors [||, || 
studying the mean field version of the Ising spin glass, the Sherrington-Kirkpatrick (SK) model 
jjj, without using replicas. After the solution of the SK model by Parisi ||, the marginality 
of one of the most dangerous eigenvalues across the whole spin glass phase was proven by two 
different technics |9, 1C[. In the truncated version of the SK model, valid close to T c , a whole band 
of zero and near-zero eigenvalues were found [11]. The nature of the spin glass phase in finite 
dimensional models with short range interaction has been the subject of intensive debate, and 
the question has not been settled until now. Nevertheless, both conflicting theories, "droplet" 
picture or phenomenological scaling on one side [12, OL [L4| and the ultrametrically organised 



complex phase space structure of Parisi on the other pi fjj^ 16 1, predict a marginal spin glass 
phase with an infinite spin glass susceptibility. 

To give another example of a quenched disordered system with a massless low temperature 
phase, we can mention the long-range correlated random manifold problem in d = D + N 
dimensions, where D is the intrinsic dimension of the manifold. For ./V — > oo, while D < 4 is 
kept fixed, the zero momentum limit of a family of eigenvalues of the mass operator goes to zero 
[17|. (In the notation of Ref. [17], \f = o(x;x,x) = 0, < x < 1, where the momentum vector p 
is D dimensional.) 

The idea that a continuous symmetry emerges in the Ising spin glass, which is a discrete 
model, giving rise to Goldstone modes was raised twenty years ago [18]. Bray and Moore argued, 
that it is the replica limit n — » 0, which may be responsible for this rather strange phenomenon. 
In Sec. ^, we will recall how the permutation invariance of the n replicas becomes a continuous 
symmetry in systems whose ordered phase is organised in the hierarchical way proposed by 
Parisi, but only in the limit R — > oo, where R is the number of steps in Parisi's construction. 
All the examples above belong to this class, except the "droplet" theory, which corresponds 
to a replica symmetric (R = 0) picture. The symmetry argument we will use applies for a 
replicated system, and is a direct consequence of the replica trick. How to derive Goldstone's 
theorem in the original models and what is the continuous symmetry there, is not clear for us 
and remains an open question. The results of Sec. ^ were derived earlier using different, though 
closely related, infinitesimal permutations Jl9[ . 

The main result of this paper is left to Sec. |3|, where it is shown how the somewhat misterious 
reparametrization invariance p5j follows from the permutation symmetry of the n replicas. Ward 
identities can then be derived, at least in principle, from the disappearance of the second, third, 
etc. variations of the free energy funcional. We have thus two different methods to obtain Ward 
identities: invariance under infinitesimal permutations and/or reparametrization. To make the 
presentation clear, computational details are left to the Appendix. 
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2 Permutation symmetry and Ward identities 



Our starting point is the free energy F, expressed as a functional of the order parameters q a p, 
a, f3 = 1, 2, . . . , n and a < (3, with n the replica number. The order parameters q a p are usually 
regarded as the elements of a real, symmetric matrix, zero along the diagonal (q a p = qp a , 
Qaa = 0)j f° r our purposes it will be useful to think of them as components of a ^n(n — 1) 
dimensional vector q. The association between the matrix q a p and the vector q is obvious: one 
lists the elements of q a p with a < (3 in any prescribed order (say row by row) to form a column 
vector. 

Although, the symmetry arguments we are going to use are completely general, we wish to 
give here some specific examples. The reader may think of F{q) as 

• the functional appearing in the integral representation of the quenched averaged free energy 



of the SK model |16|; 



• the Legendre-transformed free energy, with respect to a source h a g , of a replica field theory 
with a Lagrangean C((p a s) invariant under any permutation of the n replicas; 

• a functional obtained by a second Legendre-transformation [^] of a replica field theory 
with a Lagrangean £(4> a ) (the source h a p now couples to products like 4> a (fip). In this case 
the notation G a p is preferred to q a p. The random field Ising model [^] and the random 



manifold problem [17, [2^] are good examples where such a functional has an important 
role. Although the diagonal elements G aa are no longer zero now, this leads to only slight 
modifications in the arguments and has no influence on the results. 

By construction, the free energy is invariant with respect to the permutations P of the 
replicas: 

F(q')=F(q) for q' a/3 = q P<xPfj . 

Permuting the components of a vector will not change its length, so P generates an orthogonal 
transformation: 

q' = Oq , T = 1 . 

The invariance of F under O implies that its gradient transforms as a vector: 

Throughout this paper, Eq. || will be used for a vector q a p built up by Parisi's hierarchical 
construction || 16 1, and having the following properties: The (for the time being) discrete values 



of the matrix elements will be called q r , and the sizes of the hierarchically arranged blocks p r , 
r = 0, 1, . . . , R+l and, by convention, pq = n and pr+i = 1. r is the overlap of the replicas a and 
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j3, r = an/?, i.e. q a p = q r , and by definition atla = R+l. Both series of parameters are assumed 
to be monotonic: q s < q r and p s < p r for s < r. R is the number of replica symmetry breaking 
steps; our main concern will be the evolution of the symmetries of the system as R — > oo. For 
large but finite R the q's and the p's fill the intervals [qo,qii] with qo > 0, qn < 1, and [pi,Pi?J 
with pi > po = n > 0, pr < 1, respectively, in a quasi continuous manner: 

g r+ i-g r = 0(l/#) , r = 0,l,...,R-l , 
q r +i + q r -l-2q r = 0{1/R 2 ) , r = 1, 2, . . . , R - 1 , 
p r+ l-p r = 0(1/ R) , r= l,2,...,i?-l , 

Pr+l + Pr-l ~ 2 Pr = 0(1/ R 2 ) , T = 2, 3, . . . , R - 1 (2) 

(the difference pi — po and p_r+i — pr may be of O(l)). A vector q associated with a matrix q a p 
with the above structure and continuity properties as in Eq. ^ will, in the following, be called a 
Parisi vector. Note that we have not assumed the g's and p's to be stationary. In fact, in most 
of what follows we will be considering symmetries that are present irrespective of whether we 
are at a stationary point or not. 

Let us consider now the action of a special permutation of replicas on a Parisi vector. The 
permutation will be chosen in such a way as to interchange two blocks of size p r +2 and leave 
the rest unchanged. The replicas belonging to these two blocks will be labelled by a. L and 
respectively, i = 1, 2, . . . ,p r +2- The permutation in question will then act as P ai = and 
Pfc = on for i = 1, 2, . . . ,p r +2, and as P a = a for a outside the two selected blocks. If the two 
blocks belong to the same block of size p r +i, i.e. if aj fl = r + 1, this permutation is just an 
element of the residual symmetry group that remains after replica symmetry breaking (RSB), 
and will leave the Parisi vector invariant (q' = q). If the two blocks are chosen farther apart 
(«i H pi < r + l), however, q' will not be a Parisi vector any more, and, depending on P, a 
smaller or larger difference q' — q develops. The smallest change is expected for on Pi Pi = r. For 
this choice, the only nonzero components of q' — q are the following: 

(?' - q)a a = {q - q)/3iS = q r +i - q r , 

for 

ain^f = (3in5 = r and ain5 = Pi(lj = r + l , 

and 

(?' - q)a a = - <?W = q r - q r +i , 

for 

airi7 = /3jn<5 = r + l and a;ri(5 = ftn7 = r . 
The length of q ' — q is 



\J{q ' ~ q} 2 = J^2(Q' = \j^Pr+2{Pr+l -Pr+2)(q r +l ~ 9r) 2 (3) 
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which is of 0(1/ R 3 ^ 2 ), i.e. infinitesimal. 

The idea of constructing an "infinitesimal" permutation has been around for a long time in 



the replica approach to random systems. The first such transformation appeared in [23|, where 
it was shown that a suitably chosen linear combination of replicon eigenvectors added to a Parisi 
vector will result in a "reparametrization" of this vector, i.e. in a modified sequence of param- 
eters q r and p r . This observation will be further developed below. Shortly thereafter Goltsev 
introduced a set of "infinitesimal permutations" and showed how any finite permutation can be 
built up from infinitesimal ones |24|], without exploiting the full potential of these infinitesimal 
generators. Independently, Parisi and Slanina rediscovered the same transformation in a random 
polymer context p5[| . Finally an infinitesimal transformation closely related to the one above 
was used by ourselves to derive Ward identities for the spin glass [jnj, in the rest of this section 
we recall some of the results derived thereEI. 

After finding infinitesimal symmetry transformations, we can follow the usual steps to obtain 
Ward identities in a system with a continuous symmetry ||. Since q' is very close to q we can 
expand the left hand side of Eq. |] to get 

Of-f = M(q'-q) + ... (4) 

where / is the gradient vector 

- dF 

and M is the mass operator with components 

8 2 F 

M »Pr/S = o o • (6) 
oqapdq^s 

Evaluating the derivatives in Eqs. || and || at q a p (a Parisi vector) we get a gradient vector / with 
the same Parisi-like structure, and a mass operator which has the structure of an ultrametric 
matrix. (Ultrametric matrices commute, by definition, with all the elements of the residual 
group. Their structure was studied in detail in pq| .) Eq. || can then be analysed by the block 
diagonalization procedure described in |2(|. As shown in that paper, the relevant quantities 
representing the diagonal blocks are: 

• the one-dimensional replicon (R) blocks (i.e. eigenvalues) A(r; k, I), r = 0, 1, . . . , R; k, I = 
r+l,r + 2,...,R+l, and 

• the (R + l) x (R+ l)-dimensional longitudinal-anomalous (LA) "kernels" of the LA blocks, 
Kf,(r, s), r, s = 0, 1, . . . , R and k = 0, 1, . . . , R + 1. 



1 In Ref. |19|, at editin g, the equation nu mbering was messed up. The puzzled reader will find the correct 



numbering in the original :ond-mat/9802166 version 
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The explicit expressions for X(r;k,l) and K(r,s) are given in p6f| . After block diagonalization 
the various components of Eq. |4| give, to leading order in 1/R, the following set of equations: 

fr+i ~ fr = (q r+ i ~ ? r )A(r + 1; r + 2, r + 2) + 0(l/i? 2 ) , (7) 



A(r + l;r + 2,r + 2)-A(r;r + l,r + 2) = 0(1/ R) , 
A(r;r + 2,r + 2) - A(r;r + l,r + 2) = 0(l/i?) , 

K r+2 (r,s) -if r+2 (r + l,s) = 0(l/i?) , (8) 

with r = 0, 1, . . . , i? — 1 and s = 0, 1, . . . , R. Eq. ^ expresses a continuity property of the mass 
operator, while Eq. |?| is a Ward identity which establishes a relationship between the first and 
second derivatives of F. In the limit R — > oo, after introducing the continuous parameter x by 

* = ^Y 0<x<l , (9) 

(0) becomes 

df(x) dq(x) 

— ; — = — ; — Mx;x,x) , < X < 1 . (10) 

ax ax 



If, finally, we take Eq. 1C at a stationary point where / = 0, we obtain 

X(x;x,x)=0 (11) 

for all x where dq ^ ^ 0. (If, as is often the case, there is a breakpoint x\ beyond which q(x) 
and also \(x;x,x) are constant, ( |TTj ) still holds in the limit x — ► x^[ , and, by continuity, it also 



holds for x > x\.) The status of Eq. |ll| is that of a Goldstone theorem for spin glasses: under 



the assumptions of permutation invariance and the existence of the continuous limit R — > oo, 



(11) follows, independently of the concrete form of the free energy functional F. 

It is straightforward, at least in principle, to apply the above ideas for deriving higher order 
Ward identities. The mass operator M is the second derivative of the free energy, and thus it 
transforms as a tensor under an orthogonal transformation O: 

M' = OMO T , (12) 

where M' refers to the derivatives in (p) evaluated at q' = Oq. Expanding the left hand side of 
Eq. [l^ around q, and introducing the 3-point vertex function as 

d 3 F 
dq a pdq jS dq^v 
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we get 

(OMO T ) aW - (M) Q/3)7 5 = Wap, y s,p, - + ■ ■ ■ ■ (14) 



Eq. [L4| is the starting point for obtaining the set of Ward-identities establishing the exact connec- 
tions between components of the 2- and 3-point vertices. To accomplish this work, a procedure, 
similar to that of block diagonalisation of the 2-point vertices, transforming the 3-point vertices 
into a "canonical" form is needed. A simple example of such a relationship between 2 and 



3-point vertices is worked out in [19]. 



3 Reparametrization invariance -a subtle corollary of permu- 
tation symmetry 

As we see now, the symmetry responsible for the Goldstone modes is permutation invariance 
that becomes "continuous" in the limit R — > oo. Nevertheless, in Ref. [la] two of us purported to 
derive (|ll]) from a rather liberal interpretation of an other kind of symmetry: reparametrization 
invariance. The obviously correct result was obtained by using some uncontrollable calcula- 
tional steps, such as forgetting about the differences between derivatives with respect to x or 
p{x). Having two independent symmetries leading to the same exact property of the system, 
i.e. the masslessness of some of the modes, is rather unlikely, and one can suspect that a con- 
nection between them must exist. In this section, we will find this connection showing how 
reparametrization invariance follows, in the limit n — > and R — > oo, from the primary symme- 
try of the equivalence of the n replicas. 

The first thing we need is a clear definition of the reparametrization transformation. It 
connects two Parisi vectors q —* q' with the following properties: 

Pr = Pr+Spr , T = 1, . . . , R \ (15) 

q' r = Q(Pr) , r = l,...,R ; (16) 

q'o = <?o , (17) 

supplemented by the condition that the endpoints be invariant, 5p\ = 5pn = 0. The function 
Q(p) is obtained from q(x) and p(x) by eliminating x between them, and it remains fixed under 
a chain of reparametrization transformations (see Fig. ||) . 

At this point we wish to stress that the reparametrization invariance is not a property of the 
free energy functional F(q a p), but of the functional F(q r ,p r ) obtained from F by restricting its 
argument q a p to the subspace spanned by the Parisi-like vectors. This subspace may be called 
longitudinal, and, by extension, F may be called the "longitudinal" free energy. F remains 
unchanged when moving from q to q' by the reparametrization transformation of Eqs. |l5|]l7|: 



F(q' r ,p' r ) = F(q r + 5q r ,p r + 5p r ) = F(q r ,p r ) . (18) 
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Q(P) 




1 1 1 1 

p. P R 1 p 

Figure 1: Q(p) defines a chain of reparametrization transformations 
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In the above equation q and q' are Parisi vectors with the properties given in Eq. ^, it expresses, 
therefore, an exact symmetry of F only in the limit R — > oo. For finite R the invariance of F 
is only approximate, valid within corrections of 0(1/R). Eq. 18 can be cast into an equivalent 
form using variations of F: 

5 {k) F[q(x),p(x)} = (19) 

for any positive integer k. In Eq. [19] only p(x) is varied independently, with the constraint of 
fixed endpoints 5p(0) = 8p(l) = 0. By Eq. [l6|, 5q(x) can be expressed in terms of 5p(x), an 
expanded form of which is useful to be recorded here for further reference: 



6q(x) = Q'\p(x)]6p(x) + -Q"[p{x)]5p{x) 2 + ... 



(20) 



To prove the invariance properties in Eqs. [T^ or [l9|, AF = F(q' r ,p' r ) — F(q r ,p r ) will be 
expanded up to second order in the 5q r J s and <5p r 's. A simple Taylor expansion yields 



AF 




d 2 F d 2 F d 2 F 

-5q r Sq s + 2- — - — 5q r 5p s + - — - — 5p r 5p s 



dq r dq s 



dq r dp s 



dp r dp s 



(21) 

On the other hand, the displacement between the two "longitudinal" (Parisi- type) vectors in the 
full vector space of the (7 Q/ g's gives an equivalent form for AF. Using the short-hand notation of 
brackets as in the Appendix (Eqs. |33| , |34| ), we can write 



AF=(f\q' -q) + -(q' -q\M\q' -q) + . 



(22) 



/ and M were defined in Eqs. [5] and |(| and under any permutation P of the replicas, corre- 
sponding to an orthogonal transformation O in the order parameter space of the <? Qij g's, they 
transform according to Eqs. |l| and 12, respectively. As pointed out in Sec. ^, this is a clear 
consequence of the permutation invariance of the free energy F{q a p). Furthermore, taking the 
permutation P from the residual symmetry group corresponding to the ultrametric construction 
of the vector q (which is obviously different from the group of transformations defined by q'), 
i.e. Oq = q, we get 

/' = / = Of and M'O = MO = OM. (23) 
Thus / and M are ultrametric, and Eqs. 



|, p9| of the Appendix can be used in Eq. |22|. Putting 
everything together from the Appendix (especially from Eqs. EH1 32, 38 and HO), a comparison 



of Eqs. |2J and |22j yields: 

'-{p r -p r+ i)f 
'-(q r - q r -i) 



dF 
dq r 
dF 

dp r 



n 
2' 
n 
2 1 



fr ~ \{lr ~ q r -i)\(r;r + l,r + 1) 



0,...,R 
1,...,R 



(24) 
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From this, in the limit R — ► oo, the first variation 8^'F follows immediately. Introducing the 
continuous variable x (see Eq. 0): 



S^F = - dx [q'(x)f{x) - Q'\p{x)]p'{x)f{x)] Sp(x) 
2 Jo 



(25) 



in deriving the above expression, use has been made of Eqs. 21 and 24, together with the 
reparametrization constraint Eq. p(1. 



The integrand in Eq. 25 is identically zero, which is a direct consequence of the definition 
of Q(p): Q[p(x)] = q(x). Thus the disappearence of the first variation, 5^F = 0, seems to 
be somewhat trivial. We must emphasize, however, that it is true along the whole chain of 
reparametrization transformations; a result which follows from the crucial step in the proof: 
if F(q a p) is invariant for any permutations of the replicas, then its derivatives taken at an 
ultrametrically structured q are themselves ultrametric. The permutation group of the n replicas 
is "large" enough to include all the subgroups defined by the Parisi-type vectors along a path 
of consecutive infinitesimal reparametrization transformations. 8^'F is, therefore, identically 
equal to zero, implying the two equivalent form of reparametrization invariance, Eqs. [l^ and |i~9|. 

With its proof now accomplished, we can proceed and use Eq. ^ for k = 2, 3, . . . to derive 
Ward-identities (in principle to any desired order). They must not be different from those fol- 
lowing from the primary symmetry of permutation invariance using infinitesimal permutations, 
as in Sec. |2[ To show this, we compute 5@'F. The second partial derivatives of F can be 



calculated, by intensive use of the Appendix, just like the first ones were in Eq. 24. The results 
are as follows: 



d 2 F 

dq T dq s 

d 2 F 
dp T dq s 



d 2 F 
dp r dp s 



Tl ft 

-(p r -p r +i)(Ps -p s +i)K (r,s) + -(p r - p r+ i)X(r; r + l,r + 1)4 
r,s = 0, ...,R , 

j(q r ~ Qr-i)(Ps ~ p s+ i)K (r, s) + 6% + 



Kr 



n 



[-fr + (Qr ~ g r -i)A(r; r + 1, r + 1)] 



(26) 



r = l,...,R and s = 0, ...,R 



n 



-(q r - q r -i)(q 8 ~ q s -i)K (r, s) 
r, s = 1, . . . , R . 



It is now straightforward to compute 5( 2 'F using Eqs. 2C ,21,24 and ^6|. Note that the second term 



in Eq. 20 multiplied by the first partial derivative of F with respect to q(x) also contributes. 
The rather lengthy calculation consists of mainly partial integrations, surface terms always 
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disappearing because of the constraint 5p(0) = 5p{l) = 0. After a lot of simplification, the 
Ko(x,y) terms cancel each other, yielding the surprisingly simple result: 



6&>F = £ f 1 dx [q'(x)X(x; x, x) - f(x)} Q'[p{x)} 5p{xf . (27) 
£ Jo 

Since 5p{x) is arbitrary, and assuming a monotonic Q(p), the vanishing of the second variation 
is equivalent with the vanishing of the expression in the brackets in Eq. |2?]; i.e. Eq. 10 is now 
regained from the invariance of the free energy F when the "gauge" p(x) is changing. 

An example of two equivalent "gauges" with fixed endpoints is shown in Fig. |!| We can even 
deform curve (b) to curve (c) in Fig. |3| with the plateaux around x = and x = 1. Supplementing 
the definition of Q(p) by the plateaux regions drawn by dashed lines in Fig. [j], we can easily 
figure out that the parametrization of curve (d) represents the same Parisi vector q a p as curve 
(c). In that case, q(x) has also two plateaux with <5 m i n and (Jmax- Thus, the "longitudinal" free 
energy remains invariant for reparametrizations where the endpoints p(0) and p{l) move apart: 
< p(0) < pi and pr < p(l) < 1. The most common gauge, p(x) = x, introduced by Parisi ||, 
is also displayed as curve (e) in Fig. [|. 

As a last remark, we want to stress that reparametrization invariance expresses simply the 
fact that it is not the gauge p{x) but Q{p) which has a physical meaning. (For the Ising spin 
glass, the probability distribution of pures states is the derivative of the inverse function p{Q) ) 
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Appendix 

In the starting phase of Parisi's ultrametric construction (see the references |l6|) the block sizes 
PrS are large integers with the property that each p r +i is a divisor of p r for all r = 0, 1, . . . , R 
and the continuum of the p r 's with the inverted monotonicity of a Parisi vector, defined in Sec. ^, 
is obtained after taking n — > 0, R — > oo. (For a more mathematical treatment, see f27j.) We 
adopt a similar construction in the reparametrization procedure and consider p r +i a common 
divisor of p' r and p r , while p' r itself is a divisor of p r ; of course, this is valid at the stage before the 



"upside down" continuation of n from a large integer to zero, and 5p of Eq. 15 can be considered 
an arbitrary function of x with the only condition 5p(0) = Sp(l) = 0. 

In this appendix, however, we remain in the domain of large integer block sizes p r and p' r , for 
all r = 0,1, ... ,R and finite R, and the properties of the previous paragraph will be assumed. 
This enables us to expand the vector q ' — q, the displacement under the reparametrization 
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Figure 2: Reparametrization with fixed endpoints 
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Figure 3: Extending reparametrization transformation to gauges with moving endpoints (see 
text) 
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transformation, in terms of basis vectors. These basis vectors form a complete set in the IiilL_L)- 
dimensional space of the q a ps, and they block-diagonalize any generic ultrametric matrix. (For 
a detailed analysis of the structure of this non-orthogonal basis, see [^6|.) Following |23|], we 
can figure out that only two types of basis vectors are involved in the expansion: a replicon one 
j£{R,r) anc j a longitudinal one q( L ' r \ g( R ' r ) is just the linear combination of the (r;r + l,r + 1) 
modes which shows up in Eq. (7) of [23|, and as such, it is an eigenvector of any ultrametric 
matrix with the eigenvalue A(r;r + l,r + 1). It has nonzero components only for a Pi (3 = r, 
taking two different values depending on the overlap of a and j3 with respect to the ultrametric 
structure defined by the new block sizes p' r -s: 



q a g = 2 i a ifi m the same p' r block, 



,(R,r) (Pr ~ 2p r+1 )(p r - p' r ) 

' Q/3 2pl +1 

gg' r) = JPr- 2 Pr+lWr-Pr+l) ^ a> /Hn different ^ blocks. ( 28 ) 

As for the longitudinal vector q ( L ' r ), it has zero elements everywhere, except for a n = r: 

q% T) = 1, «n/3 = r; 

q% r) = 0, an(3^r. (29) 

r = 0, 1, . . . , R in Eq. |2^, while, remembering that p' = po = n, the replicon vector of Eq. [2^ is 
not defined for r = 0. 

The displacement vector q' — q can be expanded as 

q' -?=E (^ (R ' r) ^ (R ' r) + ^ (L ' r) ? (L ' r) ) + (flfi " ») ? (L ' 0) > (3°) 

with the coefficients if( R ' r ) and K( L > r ) determined from the conditions 

(q f — q) a f3 = q' r — q r ck, /3 in the same p' r block, 

(<7 ' ~ <T)af3 = q'r-i — Qr a, in different p' r blocks, 

for r = a fl = 1, . . . , R. (The r = case is trivial, with (<f' — q) a p = q'o — 9o> leading to the 
simple last term in Eq. [!(].) Using Eqs. 28,29, it is straightforward to obtain i^( R > r ) and K^ L ' r \ 
For later reference, it is useful to express them in terms of 5q r = q' r — q r and 5p r = p' r — p r : 

= -, r [(q r ~ q r -i) + Sq r - 5? P _i] , (31) 

(Pr ~ 2p r+ i ) (pr - p r+ i ) 

j£( L ' r ) = fiq r — SlJ. fip r (5q r -\5p r — 5q r 5p r ) . (32) 

Pr ~ Pr+1 Pr ~ Pr+1 
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What we need in the main text, is the scalar product 



{f\q'-q) = ^2 /a/ 3 - 9)0/3 



and the matrix elements 

(<?' - g I M I g' 



E ~ 9) a/3 M aPa5 (q' - q) 7S 

a<f3,-y<8 



(33) 



(34) 



for an ultrametrically structured vector / and matrix M. Since / is now a longitudinal vector 
(orthogonal to any replicon one), and also using Eq. it follows: 



</|<? (R ' r) ) = , 



(/k (L ' r) > = S(Pr-Pr+l)/, 



Combining Eqs. 3C and 35, 



(f\q'-q) = ^ £(Pr " Pr + l)frK^ r) + £(n - Pl )f 8q 



(35) 



(36) 



r=l 



obtains. To find a similar formula for the matrix element, we must use results from Ref. [|26| 



Mg (R ' r) = A(r;r + l,r + l)qr' 







(37) 



s=0 



where the block matrix element Mg°r can be expressed by the longitudinal kernel Kq(s, r) using 
Eq. 44 of Ref. M: 



Mj® = A(r; r + 1, r + 1) + ^(p r - p r+1 ) K (s, r) 



(38) 



From Eqs. 3C and 37, it is now straightforward to find the following expression for the matrix 
elements: 

(q'-q\M\q'-q) = 

Ft R 

J2 # (R ' r)2 A(r; r + 1, r + 1) (q^ \ q^} + £ i^ L ' r >i^') M$ <g( L ' r ) | gCM) + 

r=l r,r'=l 

+ <5g E ^ (L ' r) (m$ (gCM I g^)) + i<5 (^ L '°) I g( L '°))) + 5g 2 M$ <g< L '°) | q^). (39) 



r=l 
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The scalar products occuring in the above formula can be easily computed from the definitions 
in Eqs. and [29|: 
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